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Abstract. The origin of cosmic rays have remained a mistery for more than a century. JEM-EUSO is a pio¬ 
neer space-based telescope that will be located at the International Space Station (ISS) and its aim is to detect 
Ultra High Energy Cosmic Rays (UHECR) and Extremely High Energy Cosmic Rays (EHECR) by observing 
the atmosphere. Unlike ground-based telescopes, JEM-EUSO will observe from upwards, and therefore, for 
a properly UHECR reconstruction under cloudy conditions, a key element of JEM-EUSO is an Atmospheric 
Monitoring System (AMS). This AMS consists of a space qualified bi-spectral Infrared Camera, that will pro¬ 
vide the cloud coverage and cloud top height in the JEM-EUSO Eield of View (EoV) and a LIDAR, that will 
measure the atmospheric optical depth in the direction it has been shot. In this paper we will explain the effects 
of clouds for the determination of the UHECR arrival direction. Moreover, since the cloud top height retrieval 
is crucial to analyze the UHECR and EHECR events under cloudy conditions, the retrieval algorithm that ful¬ 
fills the technical requierements of the Infrared Camera of JEM-EUSO to reconstruct the cloud top height is 
presently reported. 


1 Introduction 

Although cosmic rays have been discovered over one 
century ago, still their origin remains unknown. JEM- 
EUSO (Extreme Universe Space Observatory on-board 
the Japanese Experiment Module) is a pioneer space- 
based experiment meant to detect Ultra High Energy Cos¬ 
mic Rays (UHECR) and Extremely High Energy Cosmic 
Rays (EHECR). It will observe the Extensive Air Showers 
(EAS) produced in the atmosphere when the UHECR and 
EHECR interact with the atmospheric nuclei. 

To estimate the physical properties of the UHECR 
such as energy, arrival direction and its composition, the 
EAS fluorescence light profile and the EAS Cherenkov 
light profile has to be measured. Since JEM-EUSO will 
cover a very large observation area in the atmosphere 
(~ 10^^m^), a properly monitoring of the atmospheric con¬ 
ditions is mandatory. The Atmospheric Monitoring Sys¬ 
tem (AMS) Q of JEM-EUSO will consist of a bi-spectral 
and space qualified Infrared Camera and a LIDAR (Light 
Detection And Ranging). The bi-spectral IR camera gives 
the cloud coverage in the EoV of the main telescope. With 
the difference in the brightness temperature in the IR bands 
and a radiative retrieval model, the clouds altitude is pro¬ 
vided. Thanks to the laser back-scatter signal, the optical 
depth (r) profiles of the atmosphere in selected directions 
are measured (detecting aerosol and cloud layers). 
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To increase the statistics of the detected EAS, JEM- 
EUSO may use certain triggered events observed in cloudy 
scenarios m The main criteria to classify cloudy events 
as "reconstructible events" is the detection of the depth of 
maximum development. The shower maximum will be ob¬ 
servable in cloudy conditions if r < 1 (the cloud is opti¬ 
cally thin) or if H^ax > He (i.e., the altitude of the shower 
maximum is higher than the cloud top height) HQ) 

In this work we focus into two main features of the 
EAS reconstruction in cloudy conditions: the EAS geom¬ 
etry analysis and the reconstruction of the cloud top height. 

2 Shower geometry in cloudy conditions 

EAS signal is detected as a spot moving on the focal sur¬ 
face of the JEM-EUSO main telescope. To analyse the 
direction of the incoming shower, the direction can be un¬ 
derstood as the composition of two projections: the azimu- 
tal angle, contained in the focal plane, and the zenith an¬ 
gle, contained in the track detector plane (the plane which 
contains the shower track and the detector). The former is 
calculated with the projection of the image of the shower 
in the focal surface. The latter, from the timing informa¬ 
tion and arrival angle of the EAS photons to the detector 
0 . 

The timing fit method is based on the viewing angle 
(a), which is the angle between the direction given by two 
pixels, and the arrival time (f/). The zenith angle (6) can 
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be calculated thanks to this timing fit (Figure [T]): 



Figure 1. Scheme of the geometrical relations used for the 
shower geometry. 




where Dje is the distance from JEM-EUSO to the 
ground, Di is the distance between the beginning of the 
shower to the telescope, D 2 is the distance from the tele¬ 
scope to the EAS core, L is the shower length, y is the 
angular distance between the center of the FoV and the 
EAS core, and ti is the time at which the photons have 
been created. From the scheme [T]it is clear that: 


D2 = 


Dje 

cos(y) 


( 1 ) 


For simplicity we assume that the first EAS photon is 
emitted at =0 seconds. Then, the relation between the 
time and the zenith angle will be given by: 


(^1-^2) = 


D2 + / — Di 

c 


( 2 ) 


Which can be expressed as a function of 7 , a and 0\ 


(^ 1 -^ 2 ) 


i 

c cos(y) 



sin(a) - sm(y -\- 0) 
sm(a -\-y 0) 


( 3 ) 


For showers landing on the center of the FoV ( 7 = 0 ), 
this equation is simplified to Eq. 


D JE 

(^ 1 -^ 2 )=—- 

C 


1 -f 


sin{a) - sin{6) \ 
sin{a + 6) ) 


(4) 


In Figurej^it has been represented the viewing angle as 
a function of the arrival time for a proton induced shower 
with a zenith angle of 60° and an azimutal angle of 45° 
landing on the center of the FoV in three different cases: 
in clear atmosphere (upper panel), in presence of an op¬ 
tically thin cloud (intermediate panel), and in presence of 
an optically thick cloud whose cloud top height is slightly 
below the shower maximum (lower panel). The signal has 
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Figure 2. Viewing angle as a function of time for three differ¬ 
ent cases: in clear atmosphere (upper panel), in presence of an 
optically thin cloud (intermediate panel), and in presence of an 
optically thick cloud (lower panel). The slope of the fit is the so 
called angular velocity. The color palette indicates the number of 
photons detected by each pixel. 


been fitted with the model For the optically thin cloud 
case, even if some few photons have been scattered due to 
the cloud and their trajectory has changed, the track of the 
shower is still well defined. On the other hand, in the case 
of the optically thick cloud the signal after the cloud has 
been lost. However, a good fitting can be achieved with 
the part of the signal from above the cloud. 

In Figure the projection in the focal plane has been 
represented. To obtain the value of the azimutal angle, one 
needs just to calculate the arctangent of the slope of this 
projection. As in the previous case, for the optically thin 
cloud case (intermediate panel), even if some few photons 
have been scattered, shower track projection is still well 
defined. On the other hand, in the case of the optically 
thick cloud (lower panel) the signal after the cloud has 
been vanished. However, a good fitting can be achieved 
with projection of the signal from above the cloud. 

Thus, we can assume that angular reconstruction is 
feasible for the previously called "reconstructible" cases 
(i.e., the cloud is optically thin with r < 1 , or the alti- 
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Figure 3. projection in the focal plane of the three previous 
cases: in clear atmosphere (upper panel), in presence of an opti¬ 
cally thin cloud (intermediate panel), and in presence of an op¬ 
tically thick cloud (lower panel). With the arctangent this pro¬ 
jection’s slope, the azimutal angle is derived. The color palette 
indicates the number of photons detected by each pixel. 


tude of the shower maximum is higher than the cloud top 
height, H„ax > He). 

Nevertheless, in presence of optically thick clouds the 
angular resolution is modified. Since the signal after the 
cloud is totally lost, the HAS shower track will be short¬ 
ened, and therefore the angular resolution for an EAS in 
presence of an optically thick cloud will be equivalent to 
that of an EAS in clear atmosphere with a lower zenith an¬ 
gle. On the other hand, the apparent movement of an EAS 
in presence of an optically thin cloud will not be signifi¬ 
cantly modified, since the lenght of the shower track and 
its timing does not change. Therefore, the angular resolu¬ 
tion will be similar to that of the same EAS in presence of 
clear atmosphere. 

3 Cloud Top Height Retrieval algorithm 

We can parametrize the total extinction coefficient (e) in 
terms of the Ice Water Content (IWC) and on the fixed ice 
size distribution m Being Pground the irradiance of the 


ground at a certain temperature, and jSdoud, the irradiance 
of the cloud at another temperature (following an atmo¬ 
spheric temperature profile), we can calculate the detected 
irradiance by the IR-Camera as: 


— Pground ’ (1 Pcloud ' ^ ( 5 ) 

In Eigurej^it has been plotted, for a cloud of 500 m 
thickness at 3 km altitude, how the integrated extinction 
coefficient in the near infrared varies with the ice water 
content of the cloud. Another cloud, has been modeled 
at 5 km, and the reconstruction performance is calculated 
for both cases. In this paper, we have used an US-STD 
atmosphere temperature profile. 
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Figure 4. Example of a cloud of 500 m thickness at 3 km 
altitude (Y axis). The figure represents the integrated extinction 
coefficient in the near infrared (color scale) depending on the ice 
water content (X axis) 


3.1 LIDAR simulation 


Eor the same ice water content and the same fixed ice size 
distribution, we have simulated the extinction coefficient 
for the LIDAR wavelength, and the backscattered signal 
that JEM-EUSO would detect to understand the capabil¬ 
ity of the LIDAR to measure the cloud optical depth |[6l. 
Eigure shows an example of the simulated backscatter 
LIDAR signal, at different IWC concentrations for a cloud 
of 3 km altitude and a physical thickness of 500 m. In 
our simulation, the LIDAR resolution would be + 2GTUs, 
which corresponds to an altitude of 1500 m for a LIDAR 
signal that has been shot vertically. A more comprenhen- 
sive simulation and analysis of both LIDAR, and Infrared 
Camera is under developement. 

To model the atmosphere we have considered an hidro- 
static approximation expressed in equation 


^ ^ _p 

dz g 


( 6 ) 


Where p is the density, z is the height and g is the grav¬ 
ity constant. If we integrate last equation, and replace p 
with the ideal gas law definition p = p x R/M x T, the 
pressure can be written as: 


M-g-z 

p = poxe R T 


( 7 ) 
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Lidar Signal vs layer 



Figure 5. Representation of the simulated backscatter LIDAR 
signal, at different IWC concentrations. Notice how the signal its 
spread for several GTU, this leads that at different concentrations, 
and noise levels, the algorithm to detect the cloud from this signal 
has a resolution not lower than 2GTU. 


where p is the pressure, T is the temperature of the 
gas, g is 9.80665 m/v^, M the air molar mass is 0.02897 
kg/mol, R the constant gas is 8.3144621 J/(molxk) and po 
is the pressure at sea level that corresponds to 101325 Pa 

Q. 

Therefore, we define the atmosphere pressure under 
the adiabatic aproximation as 

pxT~^ = Constant (8) 

where k is the air heat capacity ratio, fixed as 7/5, con¬ 
sidering the atmosphere a diatomic gas of only oxigen, and 
nitrogen. Combining the pressure equation in a gas envi- 
roment, like the atmosphere, with the last equation of the 
adiabatic aproximation we obtain a temperture-height re¬ 
lation. 



Figure 6. Retrieved height using the proposed abiabatic algo¬ 
rithm of a cloud at 3 Km and 5 Km. Comparison with the recon¬ 
structed height from the brightness temperature in the B 1 = 10.8 
p m, and B2=11.5 following the same atmospheric profile used 
in the simulation is presented as well. 


As stated in the previous section, the LIDAR will pro¬ 
vide us the altitude of the cloud area seen by a certain 
pixel, where the LIDAR has been shot. Then we use this 
information, together with the temperature of the cloud re¬ 
trieved by the Infrared Camera, to calculate the adiabatic 
constant that we need in our model. To find out the altitude 
of the cloud seen by the rest of the IR camera pixels, we 
apply the adiabatic model with the obtained adiabatic con¬ 
stant and we use the temperature of the different pixels. 
In the results presented in this paper, we have simulated 
a LIDAR shoot at an IWC=100mg/m^, then calculate the 
adiabatic constant and reconstruct the rest of the cloud. 

4 Summary and conclusions 

JEM-EUSO is an space based telescope which will use 
the atmosphere to detect the Extensive Air Showers 
(EAS) produced by the Ultra High Energy Cosmic Rays 
(UHECR) when they interact with atmospheric particles. 
Since UHECR have a very low flux, it is important to in¬ 
crease the exposure of the telescope as much as possible. 
In this paper we prove how the arrival direction for events 
in presence of optically thin clouds is similar to that of 
the same event occurred in a clear atmosphere. Also, it is 
explained why it is possible to obtain the arrival direction 
of an EAS in presence of an optically thick cloud, if the 
shower maximum is visible, with a slightly worse resolu¬ 
tion than the one we would get for the same event in a clear 
atmosphere. 

To use events in cloudy conditions for our analysis, it 
is extremelly important to measure the properties of the 
clouds in the EoV of JEM-EUSO. Moreover the radiomet¬ 
ric retrieval algorithm that fulfills the technical requiere- 
ments of the IR Camera of JEM-EUSO to obtain the top 
cloud height from the temperature measured by the IR 
camera, using a single LIDAR shot for calibration is re¬ 
ported. However, as showed in this analysis, it is clear that 
the adiabatic method aproximation is far more accurate, 
than the use of the brightness temperature and the atmo¬ 
sphere profile. 
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